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Abstract. The cumulative light curves of Gamma Ray Bursts (GRBs) smooth the spiky nature of the running 
light curve. The cumulative count increases in an approximately linear way with time t for most bursts. In 19 
out of 398 GRBs with Tgo > 2 s, the cumulative light curve was found to increase with time as ~ t^ implying 
a linear increase in the running light curve . The non-linear sections last for a substantial fraction of the GRB 
duration, have a large proportion of the cumulative count and many resolved pulses that usually end with the 
highest pulse in the burst. The reverse behaviour was found in 11 GRBs where the running light curve decreased 
with time and some bursts are good mirror images of the increases. These GRBs are among the spectrally hardest 
bursts observed by BATSE. The most likely interpretation is that these effects are signatures of black holes that 
are either being spun up or down in the accretion process. In the spin up case, the increasing Kerr parameter of 
the black hole allows additional rotational and accretion energy to become available for extraction. The process is 
reversed if the black hole is spun down by magnetic field torques. The luminosity changes in GRBs are consistent 
with the predictions of the BZ process and neutrino annihilation and thus provide the link to spinning black holes. 
GRBs provide a new window for studying the general relativistic effects of Kerr black holes. 
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1. Introduction 

It is suspected that spinning black holes reside in a variety 
of astrophysical sources. Frame dragging creates a special 
region called the ergosphere in which any material or en- 
ergy must rotate in the same direction as the black hole. 
The energetic reactions near the black hole maybe respon- 
sible for relativistic jets in active galactic nuclei (AGN), 
microquasars and gamma ray bursts dKoide et al., 2002| ; 
Frail et al., 2001). In GRBs the source of the enormous 



energy in gamma rays maybe the cataclysmic formation 
of a spinning black hole involving mergers of compact ob- 
jects such as neutron star (NS) binaries or NS and black 
holes (Piran, 1999; Ruffert & Janka, 1999) and also dur 
ing or after the collapse of massive stars (MacFadyen & 
Woosley, 19991; [Paczynski, 19981; IVietri & Stella, 19981; 



Reeves et al., 2002). The central engine is hidden from 



view and only gravitational radiation and neutrinos may 
escape and reach the observer directly from the engine. 
A key feature of the internal shock model is that the ob- 
served gamma rays reflect the variability of the central 
engine and the GRB duration may be determ ined by the 
engine ( Rees fc Meszaros, 1994 ; Piran, 199E ). The cumu- 
lative output in gamma rays of a burst indirectly reflects 
the output of the central engine via a relativistic jet. The 



advantage of using the cumulative light curve is that it 
reveals the trends by smoothing the spiky nature of the 
running light curve. The cumulative light curves of most 
bursts can be approximated by a linear function of time 
and GRBs may be regarded as relaxation systems that 
continuously accumulat e energy in the reservo ir and dis- 
continuously release it ( McBreen et al., 2002a ). In a rela- 
tively small number of GRBs, the cumulative light curves 
depart from linearity in a consistent way. The selection 
and properties of these GRBs are presented here in sec- 
tions 2 and 3, discussed in section 4 and summarised in 
section 5. 



2. Analysis of the light curves of Gamma Ray 
Bursts 

A large sample of the brightest BATSE GRBs was used 
with the data combined from the four energy chan- 
nels ( Fishman fc Meegan, 1995| ). The analysis procedures 
are described in detail elsewhere (Quilligan et al., 2002; 
McBreen et al., 2001 & 2002a). The full sample consisted 
of 100 GRBs with duration Tgo < 2 s where Tgo measures 
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the burst integrated count from 5% to 95% of the total, 
319 GRBs with Tgo > 2 s and a further fainter sample of 
79 GRBs with Tgo > 100 s to include very long bursts. The 
cumulative light curves of most GRBs could be approxi- 
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mated by a linear function implying constant output over 
most of the duration of the burst (McBreen et al., 2002a). 
Two significant minorities were visually identified that are 
better described by nonlinear changes in the cumulative 
count. In category A the running count increased towards 
the tallest pulse in the burst resulting in a nonlinear in- 
crease in the cumulative profile. In category B the running 
count decreased after the tallest pulse in the burst causing 
the cumulative profile to increase at a much slower rate as 
time progressed. In category A the normalised cumulative 
profile was fit by the function 



R(ti) = I^ 



(1) 



where R(tj) is the cumulative count at time tj and Ijj^jjj 
is the minimum count at tg. In category B the function 
used was 



D(ti) = Imax - c(tmax - 



(2) 



where D(ij) is the cumulative count at time tj, tmax is the 
time at the end of the fitted section of the burst and Imax 
is the value of the count at tnrax- Imax usually has a value 
close to unity. The profiles were fit using the Levenberg- 
Marquardt nonlinear minimisation algorithm. The fitted 
sections include at least 5 pulses > 5a from the running 
profile, 20% of the GRB duration and 30% of the cumu- 
lative total. In the two categories there were 19 and II 
GRBs that satisfied the selection criteria. An increase in 
number occurs when the selection criteria are relaxed. 

The median number of pulses N in GRBs with Tgo >2 
s is only 6 ( Quilligan ct al., 2002| ) and the requirement on 
N restricts the GRBs to about half of the total. The num- 
ber of pulses is important because they may originate from 
explosions in the central engine and discriminate against 
GRBs, where the emission is not well resolved into pulses 
because of the washing out of time structure in the jet 
before the gamma-ray photosphere, and interactions with 
the external medium including the effects of early after- 
glow ( Dermer fc Mitman, 1999| ). The changes in the cu- 
mulative count presented here are quite different from the 
smooth power law decays in G RBs of the FRED (i.e. Fast 
Rise Exponential Decay) type ( Giblin ct al., 2002 ) a nd the 
smoot h decays in a sample of pulses within GRBs ( Rydc 
& Svei s 



2001) 



Most GRBs with Tgg < 2 s can also be approximated 



by a linear fit to the cumulative light curve (McBreen 
et al., ^002a|). The GRBs have a median value of N 2.5 



(McBreen et al., 2001) and no short GRBs were found to 
meet the criteria used for long GRBs. 



3. Results on the light curves of GRBs 

The 19/11 GRBs that satisfied the selection for categories 
A/B are listed in Table 1. The running and cumulative 
light curves and fits to the nonlinear sections of the bursts 
are given in Fig. 1 and 2 for three bursts in each category. 
The cumulative light curve smoothes the spiky nature of 
the running light curve. In category A the highest peak 



flux frequently occurs near the end of the nonlinear in- 
crease and subsequently the bursts stops abruptly (Fig. 
la) or continues at a much reduced rate (Fig. Ic) due to 
decrease in accretion. In category B there is a decrease 
in the running light curve that is an approximate mirror 
image of bursts in category A (Fig. 2) and the nonlinear 
slowdown in the cumulative light curve results from the 
drop in amplitude of the pulses in the running profiles. 

The GRBs in Table I have large numbers of pulses. 
The errors due to the BATSE counting statistics can be 
neglected because they are very small compared with the 
large fluctuations caused by the pulses in the running light 
curves. Equations 1 and 2 were fit to the cumulative light 
curve to obtain the best fit to the data (Figs. 1 and 2). 
The values of f3 are in the range 1.8 to 2.3 (Table 1). The 
range in /3 was estimated by varying (3 from 1.5 to 2.5 
and visually estimating the range where there was agree- 
ment between the fitted functions and the cumulative light 
curves. The values are listed in Table 1 and typically are in 
the range ±0.2. The fluctuations between the fitted curve 
and the cumulative data are due to the pulses and the time 
intervals between them in the running light curve. These 
residuals are inherent to the pulsed nature of GRBs. The 
value of c (Table 1) is given for the parabolic value (3 = 2 
for all bursts to enable comparison of the normalised cu- 
mulative profiles. The errors on the hardness ratios (Table 
1) are all less than 7% with the exception of the soft GRB 
6903 where it is 50%. 

4. Discussion 

The nonlinear changes in the light curves are usually evi- 
dent in all four BATSE channels. The achromatic nature 
of the signals seem to exclude an effect due to changes 
in the optical depth. It has been suggested that some 
of the time structure in GRB light curves could be due 



to precession of the black hole ec 



disk (Portegies Zwart et al., 1999 


Fargion, 1999; 


Romerc 


2t al., 1999; 


Thornc et al., 1986 




MacFadyen & Woosley, 


1999 


). The estimated precessional period is about 1 sec- 



uator and the accretion 



ond and is too small to account for the gradual increases 
and decreases in the GRB light curves. The nonlinearities 
in the light curves could be a subtle signature of the cen- 
tral engine that only occur in some GRBs. Most models 
of GRBs favour a hyper-accretion process into a newly 
formed black hole and in these systems the nonlinearities 
might be due to systematic changes in the accretion rate 
or to general relativistic effects of a spinning black hole. 
All the evolutionary scenarios, such as mergers of com- 
pact objects and coUapsars, create rapidly rotating black 
holes because they accrete a fraction of their mass from 
a disk. However the engine must be capable of generating 
GRBs with many pulses with progressively increasing am- 
plitudes and GRBs that appear to be their mirror images 
with harder spectra. 

When a black hole swallows matter and radiation from 
an accretion disk, its mass and angular momentum evolve. 
The evolution of the black hole was first evaluated by 
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Fig. 1. The running (dashed) and cumulative (soUd) Ught 
curves of the BATSE bursts with trigger numbers a) 3105, 
b) 3860 and c) 6963 with count per 64 ms and cumulative 
count scales on the left and right vertical axes. The inserts 
give the cumulative count (solid) and the fit of the func- 
tion (dashed) for the relevent section. The start and end 
times are listed in Table 1 and ti is shifted to zero for the 
inserts. The vertical axes in the inserts are the normalised 
cumulative count. 



Bardeen (1970) following the initial work of Lynden-Bell 
(1969|)jThe spin of the black hole, as measured by the Kerr 
parameter a = Jc/GM^ where M is the mass of the black 
hole and J its angular momentum, increases nonlinearily 
with accretion from to approach but not equal 1 for an 



extreme Kerr black hole (Thorne, 1974). The equatorial 
radius of the ergosphcrc decreases from = 2GM/c^ to 
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Fig. 2. The running (dashed) and cumulative (solid) light 
curves of the BATSE bursts with trigger numbers a) 678, 
b) 4039 and c) 6694 with the same notation as Fig. 1. 



tion is required to significantly alter the value of a which 
acquires a value a ~ 0.97 when the mass of the hole has 
been accreted. Many models of GRBs require extremely 
high accretion rates and for a typical long duration GRB 
could be about 0.1 M0S~^. There could be substantial 
changes in a during a GRB and it is necessary to know 
how the power output varies with a and time. Assuming 
a constant mass accretion rate a evolves with time as ~ 



t' 



0.5 



for 0.3 < a < 0.8 and more slowly/rapidly for a > 



0.8/a < 0.3 ( iPopham et al., 1999| ). 

There are three energy extraction processes to consider 
and precise predictions are not always available because 



of the complexity of these systems (Meszaros, 2000; Rees 



GM/c ^nd the radius of the marginally stable orbit also 1999 ). In the standard accretion process, the vi scous en 



decrea ses from rms — 6GM/c'^ to GM/c'^ and additional ergy dissipa tion is balanced by neutrino cooling. Popham 



gravitg ftional potential energy becomes available. The de- 



pendance of a o n accreted mass is giv en by Thorne (1974 ) 
and plotted in Popham et al. (1999 ). Substantial accre- 



et al. (1999 ) called it neutrino dominated accretion flow 
(NDAF) and worked out its properties including varia- 
tions in the disk mass, accretion rate, black hole mass 
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Table 1. The GRBs that satisfied the criteria for category 
A (first 19 rows) and category B (second 11 rows). The 
columns refer to BATSE trigger number, the number of 



using the Kerr metric. The integrated energy deposited by 
vv — > e^e~ increases with time as when a increases 



pulses > 5 (7, the number of those pulses included in 
the fitted region between times ti and t2 with a minimum 
value of 7 s for GRB 1440 and maximum of 140 s for 
GRB 8101, the hardness ratio (HR) which is the ratio 
of flucnce in BATSE channels 3 and 4 above 100 keV to 
that in channels 1 and 2 below 100 keV, the percentage 
of the total integrated counts (%C) in the fitted region, 
the index fi and the coefficient c x 10^"* for the fit with j3 
= 2. G RBs 3247, 6353, 6903, 8101 and 7660 are from the 
fainter sample with Tgp > 100 s. 



above an initial value of 0.5 (Figs. 21 and 22 in MacFadyen 



Burst N/N' ti,t2 



HR % C 



394 

1122 

1440 

2450 

3035 

3105 

3247 

3489 

3860 

5526 

6353 

6453 

6587 

6593 

6903 

6963 

7318 

7575 

8101 



25/12 
10/8 
15/9 
22/11 
21/12 
30/30 
25/17 



0,25 
0,11 
10,17 
45,65 
5,75 
-10,30 
80,180 



4.6 
3.2 
5.3 
3.4 
4.7 
6.5 
7.7 



13/8 
34/14 

9/5 
25/12 
28/20 
20/11 
13/5 
17/10 
12/7 
30/20 
11/9 



-6tS3- 



-6:6- 

22.4 



42 
62 
70 
49 
52 
98 
79 



1.9±0.2 
1.9±0.2 
2.1±0.2 
2.1±0.3 
2.1±0.3 
2.1±0.2 
2.1±0.2 



6.8 
49.0 
147.1 
12.0 
1.1 
6.3 
0.9 



2.0_L0.2 95:^ 



-^pa 1.0±0.3 



-6,16 
-10,50 

-5,52 
3,25 
-10,14 
-25,15 
0,17 

-3,10 
135,165 
-40,100 



4.4 
1.7 
1.7 
7.6 
4.7 
1.6 
4.1 
21.8 
12.9 
10.4 



36 
53 
46 
73 
60 
47 
52 
78 
78 
95 



2.1±0.2 
2.2±0.2 
2.0±0.3 
2.1±0.2 
2.1±0.3 
2.1±0.2 
2.0±0.2 
2.0±0.2 
2.1±0.2 
1.9±0.2 



-9t«- 

7.6 

1.5 

1.4 
15.8 
10.2 

3.3 
18.4 
48.0 

9.1 

0.5 



678 

2891 

2929 

2984 

2993 

2994 

3408 

4039 

6694 

7660 

7766 



52/44 
18/14 
44/36 
22/12 
17/11 
36/22 
44/35 
33/31 
10/9 
7/5 
22/11 



1,36 
0.5,20 
12,50 
6,20 
0,20 
6,35 
10,60 
0,40 
1,16 
25,140 
0,20 



42.6 
23.1 
13.5 
14.2 
31.1 
19.3 
6.4 
21.5 
18.3 
8.4 
48.2 



83 
72 
75 
58 
77 
91 
91 
77 
85 



1.9±0.2 
2.1±0.3 
2.2±0.3 
2.1±0.2 
2.2±0.2 
2.3±0.2 
1.9±0.2 
1.9±0.3 
2.1±0.2 
2.0±0.3 
2.0±0.2 



6.2 
22.0 
5.7 
37.2 
16.4 
6.3 
3.2 
5.1 
39.5 
0.7 
18.1 



& Woosley (1999)) and including neutrinos down to the 
last stable orbit. The calculations reveal that at higher 
accretion rates the neutrino luminosity and efficiency for 
annihilation increase. These parameters also increase with 
the Kerr parameter a. As a increases the last stable or- 
bit moves inward. The neutrino emission from the higher 
gravitational potential increases both the luminosity and 
temperature and also increases the de nsity of neutrinos be 



cause of the more compact geome try ( Asano fc Fukuyama 



2001; Balmonson fc Wilson, 199£). The model predictions 



are in good agreement with the cumulative light curves 
presented here. However the annihilati on rate scales as 
the square of the neutrino luminosity ( Meszaros, 2000|) 
and the efficiency is low if the production is too gradual. 

An alternative and more efficient process for extract- 
ing energy maybe dissipation of magnetic fields generated 
by differential rotation in the disk. This extraction process 
is largely independent of a and includes a relativistic wind 
or jet from the disk and fiares (e.g. Kluzniak fc Ruderman 



1998; Recs, 1999; Livio et al., 1999). However a large en 



ergy source is also available in the spin of the black hole 
that maybe extracted by MHD coupling in the Blandford- 
Znajek (BZ) process (Blandford fc Znajek, 1977) and the 
MHD Penrose process (Koide et al., 200S). The rotational 
energy available is 0.29 Mc'^ for a maximally rotating black 
hole. The BZ process behaves as an interaction of the black 
hole with its surrounding electromagnetic field and de- 
pends directly on the spin of the black hole. The predicted 
luminosity is L ^ 10^° a2(M/3Mo)2(B/10i^G)2 ergs s'^ 
where the magnetic field is supplied by the torus. It is in- 
teresting that L oc oc t provided B remains constant 
on the horizon. This process also predicts a cumulative 
luminosity that increases ~ t^ . Furthermore there are the 
GRBs with many pulses with progressively declining am- 
plitudes (Fig. 2 and Table 1) and running light curves that 
are reasonable mirror images of the increases in output. 
These GRBs appear to be good examples of the reverse 
process, namely, the decrease in the Kerr parameter a pro- 
duced by the slowing down of rapidly spinning black holes 
by the BZ or other processes (van Putten, 2001; Lee fc 



Kim, 2002; Li, 2002). The BZ process can account for the 



and Kerr parameter. 



The properties of these disks have 
One extrac- 



been further explored (Narayan et al., 2001) 
tion mechanism is the vv e+e^ process that taps the 
thermal energy of the disk. The annihilation of neutrinos 
close to the axis of a rotating system, where there is a 
low concentration of baryons, was proposed as a way to 



produce a clean fireball (Eichler et al., 198!;; Meszaros fc 
Rees, ] 992). A number of models have been used to esti- 



mate the ne utrino emission from NS binaries and NS and 
blac k holes (Ruffert fc Janka, 1999| ; Lee, 2000 ; Janka et al 
199£ ; I ^osswog fc Davies, 2002| ). In a more complete model 



( Popham et al., 1999 ) the neutrino luminosity and energy 
deposition in the polar regions was numerically evaluated 



behaviour of the cumulative light curves of GRBs by spin- 
up and spin-down of the black hole. 

There are significant spectral differences between the 
bursts in Table 1 with the median values of the HRs of 
5.3 and 19.3 for the two categories. The remainder of the 
GRBs with Tgo > 2 s have median values of the HR of 
4.1. The corresponding median values of the peak energies 
Epcak for available GRBs are 340 keV and 956 keV for cat- 
egories A and B ( p^^loyd-Ronning fc Ramirez- Ruiz, 2002 ). 
The GRBs in category B are among the spectrally hardest 
BATSE GRBs, with the value o f Epcak generally foll owing 
the intensity during the burst ( Preece et al., 2000 ). The 
property of the central engine that maybe responsible for 
this behaviour is the spin of the black hole. The BZ mech- 
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anism operates most efficiently for a spun up black hole 
that is threaded by a powerful magnetic field, ft appears 
that this system will generate more energetic outbursts 
with spectrally harder GRBs and observed range of pulse 
properties (e.g. Norris, 2002; ^tcBreen et al., 2002a; Lee 
et al., 2000b; Ramirez-Ruiz & Fenimore, 200C; Quilligan 



et al., |2002^ |Nakar fc Piran, 2002| ). The progenitor sys 
tems should have high angular momentum and also form 
smaller mass black holes that are spun up and down more 
easily. The high spin of the black hole and angular momen- 
tum increasing outwards can stabilise the disk against ru n- 
away radial instabilities ( Daigne fc Mochkovitch, f998 ). 



The majority of GRBs with Tgo > 2 s have cumula- 
tive light curves that increase in a linear way with t. They 
do not have substantial nonlinear behaviour that meet 
the stringent criteria adopted here. There are a range of 
possibilities to account for this steady state behaviour in- 
cluding variations in the accretion rate or magnetic field 
that mask the slow changes due to a, different progenitors 
and disk emission dominating neutrino annihilation and 
the BZ process. In the later case the majority of GRBs 
are powered by emission from the accretion disk. However 
the relative contribution from MHD processes in the disk 



and spinning black hole are a matter of debate ( Lee et al 
200080 

In a few cases there is a gap in the light curves of 
the nonlinear increases and decreases where the output 
appears to be suspended for a short interval (e.g. Fig. Ic) 
and BATSE trigger numbers 2929, 3408 and 8I0I. The 
time intervals between pulses can be accommodated in 



model s of GRBs as relaxation system s ( Ramirez- Ruiz fc 



Merlori, 2001; McBreen et al., 2002a). GRBs with Tgo < 



2 s are well described by linear increases in the cumulative 
counts and no convincing cases were found with nonlinear 
behaviour perhaps because of the short duration of the 
emission and small number of pulses. 

5. Conclusions 

The results of this study of BATSE light curve of GRBs 
are presented. The analysis revealed that the running light 
curve of a minority of bursts increase or decrease approx- 
imately linearly with time. These results are interpreted 
as possible evidence for spin-up and spin-down of black 
holes during the burst. The BZ process and neutrino anni- 
hilation are consistent with these results. The brute force 
appearance of Kerr black holes being spun up an down by 
hyperaccretion in GRBs contrasts with the more sedate 



30-15 and other sources 


Wilms et al., 2001 


Miller et al.. 


2002; 


Woo & Urry, 2002 


Koide et al., 2002 


. 
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